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Abstract 

The 1/N expansion of the two-particle irreducible (2PI) effective 
action is employed to compute universal properties at the second-order 
phase transition of an 0(iV)-symmetric N- vector model directly in 
three dimensions. At next-to-leading order the approach cures the 
spurious small- N divergence of the standard (1PI) 1/N expansion for 
a computation of the critical anomalous dimension rj(N), and leads to 
improved estimates already for moderate values of N. 
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1 Introduction 



New methods for the quantitative description of thermal equilibrium as 
well as nonequilibrium aspects of critical phenomena have a wide range 
of important applications. Topical examples are the active experimental 
searches and theoretical explorations of properties of the QCD critical point 
in the phase diagram of strongly interacting matter pQ , or the description of 
the critical dynamics of Bose-Einstein condensates in laboratory experiments 
with ultracold quantum gases [2J. On extremely different energy scales 
their quantitative understanding requires a field theoretical description of the 
equilibrium and nonequilibrium properties in the vicinity of critical points 
associated to second-order phase transitions. The latter exhibit anomalously 
large fluctuations and are characterized by universal quantities such as 
critical exponents. 

There are few nonperturbative methods in thermal equilibrium that 
can describe the large fluctuations in the vicinity of a second-order phase 
transition j^j. However, the number of methods becomes even more limited 
once nonequilibrium dynamics and the approach to thermal equilibrium is 
considered. It is important to note that mean-field-type or leading-order 
large- iV approximations are insufficient. They are known to fail to describe 
thermalization and do not properly distinguish the different universality 
classes for critical phenomena. Efforts to go beyond leading order in 
a standard 1/N expansion of the one-particle irreducible (1PI) effective 
action fail because of spurious secular terms, which grow with time and 
invalidate the approximation for nonequilibrium physics. In addition, the 
1PI 1/N expansion shows a rather poor convergence in thermal equilibrium 
for moderate values of N, and at low order important quantities such as 
the critical anomalous dimension rj(N) are known to exhibit a spurious 
divergence as iV — > [I]. 

It has recently been pointed out that a controlled description of 
nonequilibrium dynamics and thermalization of quantum fields can be based 
on a 1/N expansion of the two-particle irreducible (2PI) effective action 
beyond leading order 0E10IHlEJ[ini- For nonequilibrium dynamics the rapid 
convergence of this expansion has also been observed in classical statistical 
field theories, where comparisons with exact results are possible This 
nonperturbative approach provides a promising candidate for a uniquely 
suitable description of both the nonequilibrium as well as equilibrium physics 
in the vicinity of critical points, since it provides a controlled expansion 
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even in the presence of large fluctuations. An important step towards 
such a conclusion is, therefore, to show that the 2PI 1/N expansion indeed 
reliably describes the thermal equilibrium properties at the critical point 
of a second-order phase transition. In this work we calculate universal 
properties, employing the 2PI 1/N expansion to next-to-leading order (NLO). 
We consider the 0(iV)-symmetric scalar N- vector model in three-dimensions. 
For N = 4 this model is expected to belong to the universality class of 
the high temperature QCD phase transition in the limit of two massless 
quark flavors and for N = 1 to the QCD critical point at high baryon 
density and temperature ^3] • in the context of Bose-Einstein condensation 
in quantum gases N = 2 characterizes the relevant universality class j2] . Here 
it is important to note that the equilibrium universality classes agree in the 
relativistic and the nonrelativistic case, and we will only consider the former. 

To show the capabilities of the method we calculate the properties of 
the theory directly at the critical temperature of a second-order phase 
transition. This is notoriously difficult within perturbative approaches since 
the correlation lengths are diverging. All the universal properties at this 
point are encoded in a single critical exponent, which can be associated 
to the anomalous dimension 77. In this case the knowledge of 77 also fixes 
all universal information about the effective potential, which encodes the 
information about the critical equation of state. We show that the 2PI 1/N 
expansion cures the spurious small- N divergence of the standard (1PI) 1/N 
expansion, and leads to improved estimates already for moderate values of 
N. We finally compare the method with related approaches that have been 
employed in the literature. 



2 Model 



We consider a Euclidean field theory for a real, TV-component scalar field ip a 
(a = 1, . . . , N) with classical action 

S[ip] = J d d x (^d^ a {x)d^ a {x) + ^if a (x)ip a (x) + (cp a (x)ip a (x)) 2 ^j , 

(2.1) 

where summation over repeated indices is implied. All information about 
the equilibrium field theory can be obtained from the partition function, or 
more efficiently from an n-particle irreducible effective action , which are 
related by Legendre transforms in the presence of sources. Here we consider 
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the two-particle irreducible (2PI) representation of the effective action [T5] . 
The most general Euclidean 2PI effective action can be written as 

T[<j>, G] = S[(/>] + ^Tr In G^ 1 + ^Tr G G + T 2 [0, G] + const . (2.2) 

Diagrammatically the contribution T2[<j>, G] is given by all two-particle 
irreducible 1 graphs with propagator lines set equal to G [T3]. The classical 
inverse propagator G^^x, y;4>) = 5 2 S[4>}/S(f) a (x)5(f)i ) (y) reads 

1-1 i~ „.. x\ — ( a a i ^,2 ^ 

6iV 

A 



^O.ofcfo J/5 0) = ( + 771 + g]y c (x)0 c (x) J <5 afe <5 (x - y) 

a (x)0 fe (x)5 d (x-y). (2.3) 



3iV 

In the presence of an external source J a {x) coupling linearly to the fluctuating 
field, the equations of motion for (p and G are 



W =Ux) , ™ = o. (2.4) 
5(p a {x) 5G ab {x,y) 



3 Two-particle irreducible 1/N expansion 

In this work we consider a systematic expansion of the 2PI effective action 
r[</>, G] in the number of field components or powers of 1/N beyond leading 
order E] ■ We write 

T 2 [0; G] = T L 2 °[G] + G] + r™ LO [0; G] + . . . (3.1) 

where T^ denotes the leading order (LO) contribution which scales 
proportional to N, while is the next-to- leading order (NLO) 

contribution ~ JV°, and r^ NLO [G] ~ 1/N etc. For the 0(iV)-model 
these contributions have been derived in Ref. Ej. 2 The LO and NLO 
contributions, which we will consider here, read (cf. Fig. ^) 

1 A diagram is said to be two-particle irreducible if it does not become disconnected by 
opening two lines. 

2 Note that in Ref. E] a Minkowskian space-time is considered, whereas here a 
Euclidean metric is used. 



3 




Figure 1: Diagrammatic representation of the LO and NLO contributions to 
r 2 [0, G}. The solid lines and small dots represent the full propagator and 
the bare vertex respectively and the indices run from a — 1, . . . , N. Diagram 
(a) is Tg given by ()3.2j) . Diagrams (b) and (c) are r^ LO , where we have 
expressed this contribution as a three-loop diagram with an effective four- 
vertex containing a "chain" of bubbles. This form for r^ LO can been seen 
from flU with the aid of (j33J) and (pTHjl . 



r™[G\ 



4!iV 



G aa (x,x)G bb (x,x) 



^NLOr 



G] 



-Trln[B(G) 



A 



6N 



— / I(x,y;G)(j) a (x)G ab (x,y)(j) b (y). 



xy 



In the above equation we have defined 



A 



B(x, y;G) = 5 (x-y) + —G ab (x, y)G ab (x } y) 

bJM 

and the logarithm in Eq. (|3.3|) sums the infinite series 

A 



Trln[B(G) 
1 



QN 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



2 / \^ G ^{ x ^y) G ab{x : yyj y-^G a ^{y,x)G a ^{y,x)) +. 
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The function l(x, y; G) is denned by 

A XI" 
I(x, y; G) = —G ab (x, y)G ab (x, y) - — / l(x, z; G)G ab (z, y)G ab (z, y) , 

(3.6) 

and resums an infinite number of "chain" graphs, which can be seen by 
iteratively expanding ()3.6|) . The function I(x, y\ G) and the inverse of 
B(a;, y \ G) are related by 

B-\x, y- G) = 5 d (x -y)- I(x, y; G) , (3.7) 

which follows from convoluting Eq. ()3.4|) with B" 1 and using Eq. (j3.6|) . Note 
that B and I do not depend on 0, and r 2 [0, G] is only quadratic in at 
NLO. 



4 Universal behavior at the critical point 

In the following we concentrate on d — 3, relevant for high temperature 
quantum field theories or classical statistical models with three spacial 
dimensions. For this case the model is known to exhibit a second-order phase 
transition at a critical mass parameter m 2 = ml(N, A) for all N. In Fourier 
space we will use the notation J d 3 g/(27r) 3 = j^. 3 To show the capabilities 
of the method we calculate the properties of the theory directly for the 
critical value m 2 c . In this case the correlation lengths are diverging, which 
spoils standard perturbative approaches. All the universal properties at this 
point are encoded in a single critical exponent. 4 In the absence of external 
sources this exponent can be associated to the anomalous dimension r/. It 
characterizes the critical behavior of the propagator or two-point function: 

G(p)= M^r 2 ' (4,i) 

which is valid in the limit p 2 /A 2 — > + . Here A corresponds to an (arbitrary) 
high momentum scale which regularizes the theory. All universal properties 
are independent of A as is also shown below. 

3 If not stated otherwise most formulae are valid for general d and the restriction to 
three dimensions will only be relevant for some specific momentum integrals below. 

4 Note that deviations from the critical point are described by a further independent 
critical exponent. 
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We note that the knowledge of r/ also fixes all universal information about 
the effective potential U(<f>), which encodes the complete information about 
the critical equation of state in the presence of an external (static) source 
~ J a (cf. (12 .4|) and below). At the critical point one has 

valid for <ft — > with a nonuniversal proportionality constant that depends 
on the specific details of the model. Here the universal critical exponent 5 is 
not independent but related to r\ by the scaling relation jH] 

5 = (4.3) 

for d = 3. The presence of scaling relations such as ()4.3)1 is a consequence of 
the limited number of relevant parameters, which can regulate the diverging 
correlation lengths at the critical point. For the 0(iV)-model there is only 
one independent relevant parameter at the critical point for m 2 = m 2 c) which 
can either be associated to a nonvanishing momentum or source. As a 
consequence, one may either extract the universal behavior from (|4.1|) or 
from (|4.2|) . We emphasize that the scaling relation ()4.3|) is a robust property 
of the theory in the absence of additional relevant parameters. In contrast, 
the error involved in an approximate estimate for a difficult quantity such as 
the anomalous dimension is typically large. In the following we derive the 
relevant equations for both rj and 5, employing the 2PI 1/N expansion to 
NLO. We will then choose the simpler (former) case for an explicit solution 
of the equations in order to obtain a quantitative estimate of the universal 
behavior. 

Firstly, for a computation of the effective potential it is sufficient to 
consider a constant field expectation value. By virtue of O(N) rotations 
the most general field configuration in this case can be chosen as 



6N 

<Pa{x) = y — (f)5 al , (4.4) 
G ab (x,y) = dia,g{G\\(x-y),G ± (x-y),...,G ± (x-y)} , (4.5) 

where we have rescaled the field for later convenience. For N > 1 the 
stationarity condition ()2.4j) for the composite field G with ()3.2|) and (|3.3j) 
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translates into two coupled equations for the longitudinal and transverse 
components. In Fourier space one finds 

Gf(p) = p 2 + m 2 + 30 2 + ||A[ 3G||(q) + (i V-l)G ± (q)] 

- 20 2 I(q) - A [I(q) + 20 2 Gn(q) (1 - I(q)) 2 ] G„(p - q)} , 
Gl x (p) = p 2 + m 2 + 2 + ^|A[ G||(q) + (iV+ i)G ± ( q )] 

- A[i( q) + 2 ^G'||(q)(l-I(q)) 2 ]G ± (p-q)}. (4.6) 

Here the resummation function I is given by ()3.6|) for the configuration ()4.5|) . 
in Fourier space. We will discuss the resummation function in more detail 
below. 5 The effective potential U(4>) is determined by the effective action for 
a constant field: 

f/(0) = -^r[0,G(0)]U =const , (4.8) 

where V3 is the three-dimensional volume. Here G(<f>) = {Gii(0), G± (</>)} 
denotes the solutions of We then obtain from Q with (jT^ - lET3|) 

the relevant equation for determining the equation of state: 



dU{(j>) _ 6N<p 
86 ~ ~T~ 



m 2 + 6 2 + A / ^3G?,|(q;0) + (JV-l)G?j.(q;. 



-2I(q)G||(-q;0)}j . (4.9) 

We emphasize that the equations ()4.6|) and (|4.9|) are valid also away from the 
critical point of a second-order phase transition. In particular, the effective 
potential can be used to compute the different longitudinal and transverse 
susceptibilities in the phase with spontaneous symmetry breaking for N > 1: 

X11 ~ ^Z. > x± ~7zr- ( 4 - 10 ) 



5 The equation of motion for the field expectation value is not required since one 
considers U((f>) for all <f>. For completeness we note that the stationarity condition for 
the field J23 leads for ^ to 

A f d d q 



-m 2 



6N J (2tt) 



- d [3G ll (q) + (N-l)G ± (q)-2I(q)G ll (-q)] . (4.7) 
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Note that Goldstone's theorem is fulfilled for the 2PI 1/N expansion, which 
has been pointed out previously in Ref. [BJ. To compute the universal 
properties at the critical point, i.e. to obtain the exponent 5 from (|4.2|) . 
one has to specify in (|4.9J1 the critical value of the mass parameter m 2 = m 2 . 
The latter can be obtained for = from either equation of (|4.6|) by the 
condition G?jj _1 (p = 0) = Gj^p = 0) = 0. 

An alternative calculation of the universal properties at the critical 
point determines the anomalous dimension rj. For our current purposes it 
is sufficient to perform this calculation explicitly, which is a considerably 
simpler task. By virtue of O(N) rotations, at the critical point the most 
general field configuration in the absence of sources is given by 

G ab (p) = G(p)5 ab . (4.11) 

Equations ()4.6|) then reduce to the single expression 

G-i(p) =p 2 + m 2 + S(p) (4.12) 

where the self-energy contains a momentum-independent O(N ) part and 
both momentum-independent and momentum-dependent 0(1/N) parts, 

E(p) ^^Z G(q) -^/ q G(p - q)i(q) (4i3) 

and the chain sum is 

I(q) = 1 - (l + ±J G(q - k)G(k)) 1 . (4.14) 

The equivalence of the explicit form for the resummation function (J4.14j) with 
the implicit form of Eq. ()3.6|) for the configuration (|4.11|) can be observed 
by expansion in a geometric series. Eqs. (|4.12|) - ()4.14|) provide a closed set 
of self-consistent equations to determine the full propagator G(p). Fig. El 
shows them in diagrammatic form. 

According to (|4.1j) at the critical point the inverse propagator at zero 
momentum vanishes, i.e. G _1 (p = 0) = 0. Using this and subtracting from 
(I4.12J1 the same expression for zero momentum we can write 

G-Hp) = P 2 - ^ j (G(p - q) - G(q)) I(q) . (4.15) 
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(b) 



I 

(C) ♦^M 



Figure 2: Diagrammatic representation of self-consistent equations for the 
full propagator G(p) obtained from a 1/N expansion of the 2PI effective 
action to NLO for = 0. Diagram (a) is (|4.12|) . which expresses the 
full propagator in terms of the classical propagator Go (thin line) and the 
self-energy E(p) (shaded blob). Diagram (b) is (j4.1Hj) . which shows how 
the self-energy, the amputated one-particle-irreducible two-point function, 
can be expressed in terms of a one loop contribution containing the LO 
contribution, and a NLO contribution that involves the chain sum. Diagram 
(c), corresponding to (|4.14|) . shows how the chain sum can be expressed in 
terms of the full propagator. 

This equation is valid at the critical point and can be conveniently employed 
to extract universal properties. With the help of (j4.14|) this gives 

G- 1 (P)=P 2 + ^ |(G(p-q)-G(q))(l + ^G(q-k)G(k)) , 

(4.16) 

where we used that the integral involving the constant part of I(q) vanishes 
for p 2 /A 2 — > 0. The latter limit characterizes the momenta for which the 
universal exponent rj of the critical form for the propagator (|4.1jl is defined. 
It remains to be shown that (|4.1j) indeed solves the above equation in this 
limit. To show this, we insert (j4.1j) into (j4.16|) and obtain an equation for 
the anomalous dimension 77. The latter equation has a unique solution for 
given N as is demonstrated in the following. 

Using (J4.1j) the one-loop subintegral in Eq. (J4.16)) can be performed for 



O 



+ 



0*00 



+ -■ 
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-1 < 77 < 1/2 to give 6 



17/2-1 /, o\ ij/2-1 



6A V A 2 



k/A 

2\ ^7-1/2 

^(77), (4.17) 



with 



1 rfi - n)(r( i±2) r 

8^ 2 (r(l-2)) 2 r(l + r ? ) 



With this notation the equation ()4.1fjj) takes the form 

V\ M/2 = P 2 A /• // (p-q) 2 ^ 2 - 1 /q 2 y /2 - lN 



A 2 / A 2 3NA y q/A \V A 2 y VA 2 / 

This provides a self-consistency equation for 77 which may be solved 
numerically for given N and A/A. We have done this for a check which 
is discussed below. 

However, one can proceed further analytically to obtain rj. For the 
low momentum range of critical phenomena, p 2 /A 2 — > + , the remaining 
momentum integral is dominated by small q 2 ~ p 2 (cf. above that 77 < 1/2 
for the allowed range). In this limit we can, therefore, write 

(p-g^ v 72 ' 1 (±\ /2 - r 

a(v)n y q/A u a 2 ) Uv 





6 It can be conveniently obtained using the Feynman parametrization for non-integer 
exponents a, (3: 

3-= / dxidxo^fa;! +x 2 - 1) - 7-5 , - ; ; . 

7o 1211 ; (nA + a; a S) a+/J r(o)r(/3) 

The momentum integration can then be performed for r\ < 1/2, and the subsequent 
integration over the Feynman parameter for 77 > —1 with the result l|4.17|l . Note that we 
have evaluated the integral without regulator. One can show that subleading corrections 
obtained from keeping the momentum integration finite are irrelevant for the universal 
low momentum behavior. 
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One observes that the dependence on the coupling A dropped out completely 
from the equation. This is a manifestation of universality, which implies 
that the anomalous dimension for the three-dimensional 0(iV)-model is 
only a function of the number of field components N. After performing 
an elementary angle integration one finds with rj ^ and using the notation 
p = | p | /A and q = |q|/A: 



p 2-V 



(4.21) 

The remaining momentum integral can be performed with the help of 
hypergeometric functions. This is described in the appendix. The result 
can be written as a sum of an anomalous term ~ p 2 ~ v and a regular function 
F v (p 2 ) of momentum squared: 



p 2 ~v = p 2 + B(r]) p 2 ^ + F v (p 2 ) , (4.22) 



with 



B(r,) - 4t7(1-2t?) cos^vr) / J_\ 

W ~ (3-77)(2- V ) sin 2 (W2)iV + W J ' { ] 

The regular function can be expanded in powers of the critical momentum 
p — > as 

(l-i))(2-r,)_ 2 , „f. 1 



Since the behavior of F, n (p 2 ) for small momenta is ~ p 2 , this function along 
with the p 2 -term in f!4.22|) are subleading for a positive anomalous dimension. 
In this case, the p — > behavior is dominated by the anomalous term ~ p 2 ~ v . 
Comparing the left and the right hand side of Eq. 1)4.22)1 . one observes that 
r) has to fulfill 

B(t]) = 1. (4.25) 

We observe from this constraint with (J4.2H)) that the anomalous dimension 
indeed has to be positive for the allowed range rj < 1/2, and ()4.25)1 can be 
used to extract the critical exponent rj. In addition, we have numerically 
determined rj directly from Eq. ()4.19)) for selective values of N, which agree. 7 

7 We note that an approximate estimate for rj may also be obtained from enforcing that 
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Figure 3: The anomalous dimension rj for the three-dimensional 0(iV)-model 
as a function of N. The figure compares the results obtained from the 2PI 
1/N expansion to next-to-leading order (NLO) with the conventional 1PI 
1/N expansion to NLO. One observes that the expansion of the 2PI effective 
action cures the spurious N —> divergence of the 1PI 1/N expansion. This 
is reflected in the improved estimates of 77 for small N if compared to results 
from alternative methods [3]. In contrast, for large enough N both expansions 
converge to the same results and lim rj(N) = 0. 
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In Fig. 01 we have displayed the results for t] obtained from (J4.25|) as a 
function of iV (solid line). For comparison, we also show the corresponding 
results from the conventional (1PI) l/N expansion to NLO (dashed line) [3]. 
In contrast to the latter, one observes that the expansion of the 2PI effective 
action leads to a well-defined limit rj — > 0.5~ as N —>■ + . 8 It therefore cures 
the spurious small- N divergence of the standard (1PI) l/N expansion, which 
is reflected in the improved estimates of r\ for moderate values of N. Despite 
this qualitative and quantitative improvement, the NLO approximation 
within the 2PI l/N expansion still cannot compete with the more elaborate 
estimates from alternative methods for small N. For instance, the 2PI result 
for iV = 4 is 7] = 0.061, which is still about 35-40% too high compared to 
alternative estimates 

In the following we compare the 2PI l/N expansion with previous 
approximations based on an ansatz for the exact Schwinger- Dyson equation 
for the propagator. We note that the result from the 2PI l/N expansion at 
NLO ()4.25|) agrees with the expression obtained by Bray [T^] from his "Self- 
Consistent Screening Approximation" (SCSA). The latter approximation 
scheme corresponds also to the so-called "Bare Vertex Approximation" 
(BVA), which has been employed in the context of time evolution 
problems jjj. For BVA/SCSA one introduces an auxiliary field for a 
composite operator into the 0(A r )-model. In addition to the original 
propagator G, the effective theory then contains a propagator for the 
composite field and a two-point function mixing the original and composite 
field. The exact Schwinger-Dyson equations for the two-point functions are 
then approximated by keeping the interactions of the corresponding effective 
theory bare. As has been pointed out in Ref. [0], the BVA is not consistent 
with the l/N expansion of the 2PI effective action in the presence of a 
source term or for <fi ^ 0, since BVA sums NLO and only part of the NNLO 
contributions. For = 0, the 2PI 1 /N expansion to NLO and the BVA/SCSA 

the sublcading terms ~ p 2 in (|4.22|l cancel, which leads with H4.24[) to the condition 

(l-r?)(2-7?) = 6tt 2 r,A{ri)N . 

The solution of this equation and of (|4.25() agrees to good accuracy for N > 1, and have 
the same limit r\ — > 0.5~ as N — > + . However, only the leading contributions for p — > + 
were calculated properly since they determine the critical behavior, and this is described 
by the solutions of l|4.25|l . 

8 The limit N — » describes the universality class for the critical swelling of long 
polymer chains |17|. 
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are identical, which is the reason for the agreement observed above. We 
emphasize that the 2PI 1/N expansion can be systematically improved and 
going beyond this order requires to include the NNLO corrections described 
in Ref. 0E|. 

5 Conclusions 

The 2PI 1/N expansion represents one of the few nonperturbative methods 
that can calculate critical behavior directly in three dimensions. In 
particular, there are no improvement procedures involved such as Borel 
transformation and conformal mapping underlying results from expansions 
in 4 — e dimensions (cf. the first Ref. of [3]). In view of the improved 
behavior of the 2PI 1/N expansion as compared to the standard (1PI) 
1/N expansion, the former seems to provide a promising candidate for 
quantitative estimates. Here it is important that the 1/N expansion of 
the 2PI effective action can be systematically improved. A very interesting 
further step would take into account the NNLO corrections as described 
in Ref. Ej. But already the NLO approximation provides a valuable 
quantitative tool for studying critical phenomena. This concerns in particular 
real-time properties of quantum field theories. It should be stressed that 
the 2PI 1/N expansion presents, so far, a uniquely suitable method that 
can deal with both nonequilibrium El El El 13 UH] as well as equilibrium 
problems even in the presence of large fluctuations. Our results are very 
encouraging for an application of these methods to dynamical properties of 
critical phenomena. 
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6 Appendix 

In this appendix we discuss the analytical evaluation of the momentum 
integral in Eq. ([4.21)1 . Writing the equation as 

we decompose the integral into three terms. The first integral is 

Ri = — f dq oWl + £)" , (6.2) 
V Jo \ Pj 

which in terms of the hypergeometric function 

„ , , . T(c) f 1 , x^Hl - xY-*- 1 , , 

2FlMc > z) = r(b)r(c-b) J dx (i - xzf (0) 

becomes 

tf 1 " 1 If 1\ 
Ri = - q 9 2F1 ( 3 - 277, -77, 4 - 277, — J , (6-4) 
77 3 — 2?7 \ p / 

where we have simplified products of T-functions and relied upon the relation 

2 Fi(a, 6, c, z) = 2 Fi(b,a,c,z). (6.5) 
The hypergeometric function can be transformed further using 



2 F 1 (a,b,c,z) = ^}JZ? a \ (-z)~ a 2F1 (a, l + o- c, 1 + o - 6, - ) (G.6) 

r(6)r(c-o) v z ; 

+ lf\l ( rt H'' ^(M+6-c,l + 6-a,I 
r(a)T(c — 6) \ 2 

This yields hypergeometric functions suitable for a small p expansion, i.e. 

i\ *-2 V ^r(4 - 277) r(-3 + 77) 



2 Fi 3-277,-77,4-277,— = p 



v) ' V r(— 77) 

_ /3-2T7 77J3-277) 77 (3 - 277) 2 

+ P M -5 + -7^ x-P o P 

V 3 - 77 (2-77) 2 

(3-277) 0-~v) (2-77) 
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p 3 + 0(p 4 )^). (6.7) 



The second integral is 



R-2 



v 



77-1 pi 



V Jo 



dq q 



2— 2r? 



1 2 



P. 



n/2 



(6.8) 



which may be evaluated with the aid of the following relation where z > 1: 
1 / r(-l-a-26)r(l + 26) 







+ 



dxx a (l -xz? b = z- l - a 



r(l + a)r(l + 26) 



r(2 + a + 26) 
So we obtain 



z 



2b 



1 + a + 26) 



,,, , (6-9) 

r(-o) 

2Fx(-l - a - 26, -26, -a - 26, 1/z) . 



i?2 



f-r, / F (_ 3 + r? ) F ( 1 + ^ r( 3 _ 2?7 )r(i + v ) 



v 

1 1 



r(-2 + 2r]) 



r(4 - v ) 



(6.10) 



2 Fx(-3 + 7;, -t), -2 + r),p) 



rjp (3 — 77) 

Expanding the hypergeometric function in small p one finds 

2 F X (-3 + 77, -77, -2 + 77,p) = 1 - ^"^ p 

2 — 77 

77(3-77) , , (3-77) (1-/7) (2-77) 



■p + 

2 F 6 
The third integral is straightforward and gives 



do 2g 



1-7/ 



-2 

2-77 



p 3 + C(p 4 ) . (6.11) 



(6.12) 



Collecting the results one observes that terms containing odd powers of p 
cancel and up to O (p 4 ) one finds flj~Z2} with (jPHD and (JOIJ). 
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